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First astronomical unit scale image of the GW Ori triple system 
Direct detection of a new stellar companion. 
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ABSTRACT 

Context. Young and close multiple systems are unique laboratories to probe the initial dynamical interactions between forming stellar 
systems and their dust and gas environment. Their study is a key building block to understanding the high frequency of main-sequence 
multiple systems. However, the number of detected spectroscopic young multiple systems that allow dynamical studies is limited. GW 
Orionis is one such system. It is one of the brightest young T Tauri stars and is surrounded by a massive disk. 
Aims. Our goal is to probe the GW Orionis multiplicity at angular scales at which we can spatially resolve the orbit. 
Methods. We used the IOTA/IONIC3 interferometer to probe the environment of GW Orionis with an astronomical unit resolution in 
2003, 2004, and 2005. By measuring squared visibilities and closure phases with a good UV coverage we carry out the first image 
reconstruction of GW Ori from infrared long-baseline inter ferometry. 

Results. We obtained the first infrared image of a T Tauri multiple system with astronomical unit resolution. We show that GW Orionis 
is a triple system, resolve for the first time the previously known inner pair (separation p ~1.4 AU) and reveal a new more distant 
component (GW Ori C) with a projected separation of ~8 AU with direct evidence of motion. Furthermore, the nearly equal (2:1) 
H-band flux ratio of the inner components suggests that either GW Ori B is undergoing a preferential accretion event that increases 
its disk luminosity or that the estimate of the masses has to be revisited in favour of a more equal mass-ratio system that is seen at 
lower inclination. 

Conclusions. Accretion disk models of GW Ori will need to be completely reconsidered because of this outer companion C and the 
unexpected brightness of companion B. 

Key words. Stars:binaries: general,Stars: variables: T Tauri, accretion disks - Techniques:interferometric 

1 . Introduction GW Orionis is a single-line spectroscopic binary Classical T 

Tau ri Star with an orbit al period of 242 days located at « 400 

Stellar initial masses are controlled by the accretion of gas pc (Mathi eu et"ai1 ll99llPl The primary star's mass has been 

through disks onto the star in the protostellar phase. For a binary estimated to be 2.5 M G and the secondary's as 0.5 M with 

or multiple system this process is more complex because it po- an expected separation of ~1 AU (circular orbit of inclination 

tentially involves the presence of circumstellar an d circumbinary * 27°). Evidence of a third companion was found based on 

gas re servoirs that control the mass distribution riBate & Bonne j the analysis of radial velocity residuals (confirmed period of 

119971 lOchi et' al. 2005). Understanding the origin of such a dis- 3850 days; Latham, priv. comm.). The spectral energy distribu- 

tribution therefore requires the observation of multiple systems, tion (SED) of the unresolved complete system shows a large 

the characterisations of their masses, and of the nature of the infrared excess compared wi th a photosphere from a 5662 K 

accretion sources. G5 star (ICohen & Kuhill 1979b with a strong 10-microns silicate 
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Array config. 


UT Date 


Min. Base, (m) 


Max. Base, (m) 


Calibrators (diameters) 


A35B15C10 
A35B15C10 


2003-11-30 
2003-12-01 


16.9 


36.6 


HD42807 (0.45±0.05mas) 
HD 42807 


A28B10C0 

A28B10C0 

A28B05C10 

A28B05C10 

A28B05C10 


2004-12-12 
2004-12-14 
2004-12-15 
2004-12-16 
2004-12-20 


7.9 
10.3 


29.9 
28.6 


HD31966, HD42807 
HD31966, HD42807 
HD31966, HD42807 
HD31966, HD42807 
HD31966, HD42807 


A35B15C0 


2005-11-22 


14.3 


36.7 


HD31966 (0.43±0.06mas) 



Table 1. Epochs and array configurations of the observations. Adopted diameter used for calibrator HD 42807 (G2V, V=6.44, 
H=5.01): 0.45 + 0.05 mas and for HD31966 (G5V, V=6,7, H=5.2): 0.43 + 0.06 mas (from getcal/NexSci). Minimum and maximum 
projected baselines for each configuration are given as an indication. 



emission feature. An important near-to-mid-infrared dip in the 
SED has been interpreted as the signature of a circumbinary 
(CB) disk, whose inner part is c l earing up to 3.3 AU from the 
star system dMathieu et al.lll9911) . lArtvmowicz & Lub ow ( 1994) 
have attempted to simulate the interaction of the CB disk with 
the central system. Pho tometric time-dimming behaviour led 
IShevchenko et al.l (119981) to conclude that eclipses resulted from 
material around GW Ori B that covers the central star. 

An accurate model of the disk structure is not possible with- 
out good knowledge of the stellar components of this system. 
We present here the first long-baseline interferometer data on 
GW Ori using the IOTA interferometer and succesfully ob- 
tain the first reconstructed image of a T Tauri system with 
astronomical-unit resolution. 

2. Observations and data reduction 

IOTA was a three 45cm siderostats interferometer, located 
on Mt. Hopkins Arizona and operated by Smithsonian 
Astrophysical Observatory (ISchloerb et al.l I2006I) . The three 
telescopes (A,B,C) were movable along two perpendicular axes 
and could occupy 17 stations offering a way to synthetize a beam 
with a maximum 35 x 15 m aperture. This corresponds to a 
resolution of » 5 x 12mas in the //-band. Observations took 
place during period 2003 November - 2005 November. Four ar- 
ray configurations were used to pave the UV plane. The obser- 
vations were made in the H band using the IONIC3 instrument 
dBerger et al.ll2003l) . IONIC3 allowed each of the three baselines 
to be combined separately, which provided six interferograms. 
Fringe patterns were temporal ly encoded, detected, and locked 
with a PICNIC detector array dPedretti et al.ll2004[ 120051) . 

We used the d ata reduction method described in 
iMonnier et al.l d2004l) to extract squared visibilities (V 2 ) 
and closure phases (CP) from the interferograms. Instrumental 
V 2 and CP effects were removed by interspersing each science 
target acquisition with observations of calibrators with small 
and known diameters. Table [TJ shows a log of our observations. 
While the calibrators were chosen to match the brightness 
and spectral type of GW Ori as closely as possible to avoid 
instrument-induced biases, we estimated that a conservative 
induced CP systematic error of » 0.5° had to be added 
quadratically to the statistical error c omputed from the set of 
interferograms (IMonnier et al.l 120061) . This additionnal term 
is essentially dominated by the chromatic response of the 
beam combiner. Identically, our night-to-night calibration 
accuracy follow-up led us to adopt a conservative systematic 
V 2 calibration error of cr v i = 0.05 added quadratically to the 
statistical error on V2. Our calibrated data set resulted in the 
following statistically independent measurements: 27 V 2 and 
9 CP in 2003, 1 1 1 V 2 and 37 CP in 2004 and 6 V 2 and 3 CP 



measurements in 2005 . The calib rated data are ava ilable from 
the authors in the OI-FITS format dPauls et al.ll2005l) . 

The calibrated V 2 and CP data are presented in Fig. [TJ The 
conclusions that can be drawn from an inspection of these data 
include: 1) GW Orionis is clearly resolved at most of the spa- 
tial scales probed by IOTA; 2) The V 2 and CP curves show 
an oscillation pattern typical of a multiple component system; 
3) there is a clear non-zero closure phase signal at all epochs, 
which hints at non centro-symmetric emission. The contribution 
of the circumprimary disk's H-band emission to the V 2 curve are 
most probably negligible because their expected angular size is 
smaller than the close binary separation, which is already un- 
derresolved by our observations. We note that the CB disk in- 
ner rim is too far away to contribute significantly to the H-band 
emission. This inspection combined with the already known bi- 
narity of GW Orionis and the evidence of a putativ e companion 
to the spectroscopic system (Mat hieu et al.l dl99ll) and Latham 
priv. comm.) has led us to focus our analysis on the search for 
sources of point-like emission in the data. 

3. Modelling and results 

3.1. Visibility and closure phase modelling 

For the purpose of our fit we temporally binned the data by con- 
sidering three bins, i.e. 2003 (2003-1 1-30 to 2003-12-01), 2004 
(2004- 12- 12 to 2004-12-20) and 2005 (2005-1 1 -22). These three 
epoch bins allowed us to neglect orbital motions within each fit 
(even for the longest epoch, 2004, the eight days of observations 
only span 3% of the orbit, and therefore the orbital motion in- 
troduces negligible errors). We first fitted simultaneously V 2 and 
CP data with a binary model. To achieve this, we carried-out a 
X 2 fine grid computation in search for the minimum. We exhaus- 
tively searched the parameter space of flux ratios (/), compo- 
nent separations (p), and position angles (ff). The resulting fits 
at all epochs were poor and unacceptable (reduced \ 2 > 30), 
w hich lead us towards a triple system model as hypothezised 
bv lMathieu et alJ dl99ll) . All components were considered to be 
unresolved and the fine grid covered a range of realistic values. 
The fitting procedures defined all the above-mentioned parame- 
ters with respect to a reference source (A), i.e. not necessarily the 
primary in the usual sense of it. The fitting procedure revealed 
that epoch 2005, for which we had the least amount of data, al- 
lowed two equally plausible parameters for the tertiary compo- 
nent (C), and both solutions are included. We confirm that the 
triple model improves the reduced \ 2 considerably, finding 0.6, 
1.6, and 0.1 for 2003, 2004, and 2005 respectively. The model 
parameters and associated errors can be found in Table We 
find that the component flux ratio remains roughly constant for 
all epochs, which gives confidence in our identification of the 



J. -P. Berger et al.: First astronomical unit scale image of the GW Ori triple system 



3 



Y 



50 -5a 




Projected Baseline ( m) 



UT Time (Decimal hours) 




r 



Projected Baseline ( m) 



50 

U(m). 




4 


6 8 

UT Time (Decimal hours) 


10 










> 











Projected Baseline ( m) 



UT Time (Decimal hours) 



Fig. 1. Squared visibility as a function of the projected baseline (left) and closure phases as a function of UT time for the epochs 
2003, 2004, and 2005 (right). The inset shows the corresponding UV coverage. 



Epoch 


fslfA 


Pab (mas) 


e AB 


Self A 


Pac (mas) 


OacO 


x 1 


2003 


0.64+0.19 


3.3+0.6 


26+24 


0.21+0.02 


22.7±3.3 


247+8 


0.6 


2004 


0.57+0.05 


3.4+0.8 


189+24 


0.24+0.02 


18.6+1.3 


207+10 


1.6 


2005 


0.53±0.18 


3.4+0.6 


17+20 


0.27±0.1 


18.7+1.6 


182±9 


0.1 


2005 


(alternative outer solution) 


0.27±0.05 


12.2±2.4 


152+11 


0.1 


2004 (from MIRA) 


0.7 


3.5 


171 


0.32 


16 


217 


NA 



Table 2. Triple system parameters: flux ratio, separation, position angle, associated errors, and reduced x 2 as obtained from the 
fitting procedure and MIRA reconstruction software. 



components. We adopt the notation that the brightest component 
is referred to as "A", the other component in the close pair as 
"B", and the new outer component as "C." We continue our dis- 
cussion of these results in Sect. [4] 

3.2. Image reconstruction 

We used the MIRA software dThiebautll2008h to independently 
reconstruct an image out of the V2 and CP data. Only epoch 
2004 has sufficient data for such an attempt. This model- 
independent image (see Fig. [2j is found to support our choice of 
the triple system discussed in the previous section. Furthermore, 
even though the procedures are different, the locations of the 
three components from the image are found to be consistent with 
the best-model predictions (see Table[2]for a comparison). 



4. Discussion 

The system is resolved into three components, the spectroscopic 
pair (A & B) plus a new outer component (C). The inner pair 
shows a position reversal from 2003 to 2004 and also 2004 to 
2005. This is expected since the orbital period has been mea- 
sured to be « 242 days - furthermore, our observations hap- 
pen to coincide with times when the components are approxi- 



mately in the plane of the sky (circular orbit phases of -0.0, 0.6, 
0.0 in 2003, 2004, 2005 respectively). Unfortunately, this timing 
means that our visual orbit cannot constrain the orbital inclina- 
tion. On a positive side, this means that our measured separation 
is close to the actual semi-major axis of the orbit, 1 .35+0. 10 AU 
assuming a distance of 400pc and know ing that the eccentricity 
of the orbit is » (Mat hieu et alJ[l99l1) . This corresponds to a 
system mass of 5.6^ } "j M , somew hat higher than the previous 

estimate of 3.0 M<B jMathieu et al.lll991h . 

Perhaps the most surprising aspect of this inner pair is that 
the H-band flux ratio is only a bout a factor of ~2. Previous 
work (e.g., Mathie u et al.lll99lT> had estimated that the system 
was composed by two stars of respectively 2.5M and O.5M 
and concluded that the contribution of the secondary to the near- 
infrared excess was negligible. 

To investigate the link between the flux ratios and the masses 
we used the pre-mai n-sequence stars theoretical evolutionary 
models computed by ISiess et a"D d2000h . which were shown to 
be reliabl e enough to estimate masse s with a 10% to 20% 
precision (Iffillenbrand & White! |2004). We considered, as in 
iMathieu et alJ (119911) . that the effective temperature of the pri- 



2 We note that introducing the uncertainty on the distance (i.e 
400+3 'pc) causes an error of ss 20% in the mass estimate which is com- 
patible with this remark. 
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Fig. 2. Left: Reconstructed image of the 2004 epoch using the MIRA software. Three components are resolved (encircled dots) and 
identified according to their brightness (A, B and C respectively from brightest to dimmest in H). Blue crosses indicate the positions 
as obtained from model fitting. Centre: presentation of the data from Table [2] Square point represents the alternative position for 
epoch 2005. Right: Predicted H-band flux ratio as a function of the secondary mass with a primary mass of 3.6M Q . Shaded region 
is the flux-ratio constraint from our observations (epoch 2004). Vertical dashed lines provide the system inclination. Vertical blue 
lines are the secondary mass estimate corresponding to distance d « 414 + 7pc. 



mary was T e ff ~ 51QQK and the age w I0 6 yrs. We interpolate 
from the 10 6 yrs isochrones that this corresponds to a primary 
mass of ta 3.6M . The right most plot i n figu re [2] shows the ex- 
pected flux ratio from the ISiess et all (l2000l) computation as a 
function of secondary mass. The straight line materializes the 
observed flux ratio (i.e 0.57 + 0.05, 2004). Additionally, we use 
the mass ratio and mass function determined from radial velocity 
to com pute the requested orbital inclination as in iMathieu et alj 
(1199 11) . 

If we ta ke the most recent a nd precise determination of the 
distance by M enten et al.l d2007l) . i.e 414 + 7 pc, we can extract 
from figure |2]the following information: 

- Our observation can be explained by a binary system with 
a mass ratio close to unity at low inclination and no excess 
from circumstellar emission. Indeed, taking the uppermost 
distance estimate of 421 pc the observed flux ratio in the H 
band can be reproduced by a system seen at low inclination 
(i « 9°) with a secondary mass of » 3.1M , i.e a mass-ratio 
of the order of 0.86; 

- With decreasing distance the observed flux ratio requires 
the addition of H-band excess to the secondary and a slight 
increase in inclination. For example, if we take the low- 
est distance boundary of 407 pc, the expected secondary- 
to-primary ratio is m 0.29, which requires a secondary cir- 
cumstellar excess of roughly the same value m 0.28. The re- 
quested inclination is then w 11° 

Although this analysis is subject to important uncertainties, 
including the primary spectral type estimate, we tested that our 
conclusions hold at a qualitative level, i.e that the GW Ori inner 
system is probably seen at a much lower inclination (i.e ss 10°) 
with mo re circumsecondary H band emission than previously 
thought iMat hieu et al.\\199l\). Reconciling our ob servation with 
the observed eclipses by Shevchenko et al. ( 1998) is a puzzle we 
cannot solve. 

Another important discovery of our work is the existence of 
source C. We detected C in each epoch and found clear evidence 
for motion with respect to the inner pair AB. We compared our 
motion to the preliminary 3850 day orbit recently inferred from 
long-term radial velocity monitoring (Latham, priv. comm.) but 



did not find a solution that could simultaneously fit our motions 
and the spectroscopic orbital elements. If the 20 mas separation 
of AC is the true semi-major axis, then we would expect an or- 
bital period of ~3600 days, which is quite close to the period 
seen in the RV residuals. It is far from clear if this hierarchical 
triple is stable. Additional interferometric monitoring of the in- 
ner pair will allow the inclination to be measured directly and 
the matter will be settled quantitatively in the near future. 

5. Conclusion 

We have presented the first direct detection of the CTTS GW 
Orionis as a triple system: both the secondary and tertiary 
companion are spatially resolved. Our results show that the 
inner stars have H-band fluxes within a factor of 2 of each other, 
which suggests that the exact status of the pair should be revis- 
ited. The outer companion is clearly detectable, consistent with a 
roughly 3500 day orbital period, in line with hints from radial ve- 
locity residuals. GW Orionis offers a unique laboratory to study 
the dynamical evolution of a close, young and (perhaps) hierar- 
chical triple system and its interaction with a massive disk and 
envelope. Rapid progress in determining the full orbital elements 
of the three components is possible with current interferometers 
in the coming years. 
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